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We investigate strongly correlated phases of two-component (or pseudo-spin-1/2) Bose gases under 
rapid rotation through exact diagonalization on a torus geometry. In the case of pseudo-spin- 
independent contact interactions, we find the formation of gapped spin-singlet states at the filling 
factors v — fc/3+fc/3 (fc/3 filling for each component) with integer k. We present numerical evidences 
that the gapped state with k — 2 is well described as a non-Abelian spin-singlet (NASS) state, in 
which excitations feature non-Abelian statistics. Furthermore, we find the phase transition from the 
product of composite fermion states to the NASS state by changing the ratio of the intercomponcnt 
to intracomponent interactions. 
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Under rapid rotation, ultracold gases of bosonic atoms 
are predicted to enter a highly correlated regime, anal- 
ogous to quantum Hall systems This regime 
is reached as the number of vortices, Ny, in a Bosc- 
Einstein condensate becomes comparable with the num- 
ber of atoms, N. The relevant control parameter is the 
filling factor v = N /Ny . For scalar bosons, it is pre- 
dicted that the vortex lattice melts for v < 6 J|] and that 
a series of gapped uncondensed states appear at various 
integer and fractional v. In particular, the ground states 
at v = k/2 (with k = 1, 2, 3, ...)havc large overlaps 
with the Read-Rezayi states M, whose excitations fea- 
ture non-Abelian statistics for k > 2. Other quantum 
Hall states such as composite fermion states || have also 
been discussed. 

Given a rich diversity of strongly correlated physics 
in the scalar case, it is natural to ask what happens in 
two-component Bose gases, such as those made up of two 
hyperfine spin states of the same atoms [Q, ||. For inter- 
mediate rotation frequencies, a variety of vortex lattices 
have been shown to appear |], ^] . For a rapid rotation, 
two-component systems would offer an ideal situation in 
which to study the roles of (pseudo-)spin degrees of free- 
dom in the quantum Hall physics. As the ratio of the 
intercomponent contact interaction g^ to the intracom- 
ponent one g increases, the two spin states are expected 
to be entangled to form novel ground states. 

On the experimental front, the smallest filling factor 
achieved for rotating scalar bosons is v sa 500 |ll[] , which 
is far above the quantum Hall regime. The main difficulty 
in rotating a trapped gas faster is a fine tuning of the ro- 
tation frequency without passing through the deconfine- 
ment limit, and novel ideas for circumventing this prob- 
lem using anharmonic traps have been proposed [fl2| , [l3| . 
In order to create a rapidly rotating two-component gas, 
one can first rotate a one-component gas and then con- 
vert a portion of atoms into another component by ap- 
plying pulse lasers 0, §]. Furthermore, a very recent re- 
alization of an optically synthesized gauge field has 
opened up a new route to smaller v. With these develop- 
ments, there is now a growing prospect for realizing the 
quantum Hall regime of scalar or two-component bosons 



in the foreseeable future. 

In this Rapid Communication, we investigate strongly 
correlated phases of two-component Bose gases under 
rapid rotation by using the exact diagonalization method. 
For pseudo-spin-independent interactions = g, we 
find the formation of gapped uncondensed ground states 
at filling factors v = k/3 + k/3 (fc/3 filling for each com- 
ponent) with integer k. Precisely at these filling fac- 
tors, Ardonnc and Schoutens Jll| proposed a series of 
non-Abelian spin-singlet states (NASS) having a global 
SU (3)k symmetry; these states may be viewed as a spin- 
singlet generalization of Read-Rezayi \SU(2)k] states and 
host non-Abelian anyonic excitations Jlq, [L6| . We present 
numerical evidences that the gapped state with k = 2 is 
well described as the SU (3)2 state. We also discuss a 
phase transition that occurs at a particular value of the 
coupling ratio g-fi/g- We note that NASS states of bosons 
have also been discussed in spin-1 bosons jl7j and scalar 
bosons on optical lattices [|l8|. The two-component case 
studied here offers the simplest realistic setting which al- 
lows a detailed numerical identification of a NASS state. 

We consider a system of a Bose gas having two spin 
states (labeled by a =~f, |.) and rapidly rotating in a har- 
monic trap with the cylindrical symmetry about the z 
axis. We denote the axial and radial trap frequencies 
as u>u and respectively. In the rotating frame of 
reference, a system of neutral particles is mathemati- 
cally equivalent to that of charged particles subject to 
a uniform magnetic field. In analogy with the quan- 
tum Hall problem, we introduce the "magnetic" length 
i = y/ h/(2Mfl), where M is the particle's mass and f2 is 
the rotation frequency. 

As fi increases toward u>u, the particle density de- 
creases due to the centrifugal spreading of the gas. For 
Q w ojii , we may assume that the mean interaction en- 
ergy per particle (roughly proportional to the coupling 
constants times the particle density) is much smaller 
than single-particle energy- level spacings, hiu^±. The 
single-particle states can then be restricted to the ground 
state of the axial confinement and to the lowest Landau 
level (LLL) in the transverse directions |IJ 0. The LLL 
state with angular momentum L z = mh > in the xy 
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FIG. 1: "Charge" gap A(TV) [Eq. (g)] versus the filling factor 
per component, i^/2 = N/(2Nv), in the S'(7(2)-symmetric 
case g^i — g. The aspect ratio is set to L x /L y — 3/(2y / 3) 
and 3/(3\/3) for TVy = 6 and 9, respectively. These ratios are 
chosen to be compatible with the rectangular vortex lattice 
expected for large v M. 



plane is represented as u m (z) oc z m exp[— \z\ 2 /(A£ 2 )] with 
z = x + iy. Within this restricted subspace, we consider 
the interaction Hamiltonian consisting of intracomponent 
and intcrcomponent contact interactions: 

H iat = £ 3 Q £^-^)+-9nEE^M)' (i) 

where N a is the number of particles in the state a and 
zf s are the positions of such particles. The effective cou- 
pling constants in the two-dimensional plane are given 
by g a = a a y/8nh 3 ui\\/M and g n = a ny /8Trh 3 u}\\/M [jl9), 
where a a and a-j-j, are the s-wave scattering lengths be- 
tween like and unlike bosons, respectively. For simplicity, 
we assume g^ — g\{= g) in the following. 

To study bulk properties, it is useful to work on closed 
uniform surfaces such as a sphere and a torus, which 
have no edge. Here, we perform calculations on a peri- 
odic rectangular geometry (a torus) [Q, ^(1 of sides L x 
and L y , which contains TVy = L x L y /{2iiW) vortices for 
each component. This geometry can describe the central 
region of a trapped atomic gas where the particle den- 
sity is approximately uniform. There are TVy LLL states 
for each component. We introduce the filling factor as 
v = N/Ny, where TV = Nf + TVj_ denotes the total parti- 
cle number. We classify all states by the pseudomomen- 
tum K = (K x ,K y ) = 2Trh(m x / L x ,m y / L y ) with integer 
m x and m y . At v = p/q with p and q being coprimc, 
all energy eigenstates possess a trivial center-of-mass de- 
generacy of q, and the minimal Brillouin zone (without 
this degeneracy) consists of N 2 points, where N is the 
largest common divisor of N and Ny As in Ref. [|, 
we focus on non-negative values of K x and K y up to the 
Brillouin zone boundary, since states at (±K x ,±K y ) are 
degenerate by symmetry. Henceforth, we take the units 
in which h = 1 and £ = 1 . 

We first focus on the case of pseudo-spin-indcpcnclent 
interactions g-^i = g, where the system possesses an 



SU(2) symmetry. For this case, a mean field theory |J 
predicted that each component forms a rectangular vor- 
tex lattice with the aspect ratio of the unit cell given by 
\/3. In analogy with the scalar case Q], we expect to 
find a phase transition from the vortex lattice to gapped 
uncondensed states as the filling factor v decreases. Wc 
therefore set the aspect ratio of the torus to be com- 
patible with this lattice, such that this lattice can be 
detected if it appears. For Ny — 6, for example, we set 
L x /L y = 3/(2-^/3), with which the system can host 3x2 
units of the rectangular vortex lattice, with the longer 
side of the rectangular unit in the y direction. 

To search for uncondensed ground states, we introduce 
the "charge" gap 
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for even integer TV(> 4), in analogy with the studies 
of fermionic Hubbard models. Here, E(N^,Ni) is the 
ground-state energy for fixed particle numbers, TV-j- and 
TVj.. The results for the 5L^(2)-symmetric case = g are 
presented in Fig. |l|. We find upward spikes at v/2 = fc/3 
with integer fc, which indicate the appearance of gapped 
uncondensed states at these filling factors. 

The appearance of a gap at v = 1/3 + 1/3 can be 
naturally interpreted from the formation of the Halperin 
(221) state (2^, [23], which is an Abelian spin-singlet state. 
On a disc geometry, its wave function is written as 

* 221 = 11(4 - 4) 2 11(4 - 4) 2 11(4 - 4)- ( 3 ) 



i<j 



Here, as in Ref. [|, a tilde on the wave function indicates 
that it has to be multiplied by the usual Gaussian fac- 
tors for the xy plane. The contact interactions in Eq. (^) 
vanish for this wave function, and therefore Eq. (||) is an 
exact zero-energy ground state for arbitrary > and 
g > 0. On a torus, there appear triply-degenerate zero- 
energy ground states due to the center-of-mass degener- 
acy. Performing exact diagonalization up to TVy = 15 
and TVj- = TV4. = 5 with g^^ = g, we find a stable en- 
ergy gap of magnitude ~ 0.035g above these zero-energy 
states. 

We now investigate the origins of the gapped states 
at v = fc/3 + fc/3 with k > 2. Precisely at these fill- 
ing factors, Ardonne and Schoutens [fl5| proposed 11011- 
Abelian extensions of the Halperin (221) state, termed 
the S/7(3)fc states |24|] . On a disc, their wave functions 
are written as 
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In this construction, the TV = rk + rk bosons are first 
partitioned into k groups, each with r particles in each 
spin state f, I- For each group we write a Halperin \f 221 
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FIG. 2: Energy spectra versus the pseudomomentum K 
at the filling factor v = 2/3 + 2/3 for the S(7(2)-symmetric 
case g^i — g. The ground-state energy is subtracted from 
the spectrum. For Nv — 6 and 9, the same aspect ratios as 
in Fig. hi are used; for Nv = 12, we set L x /L y — 4/(3v / 3)- 
Upward (A) and downward (V) triangles, representing i\>\ and 
if)2 respectively, indicate the two lowest-energy states in the 
sector with K = 0, p v = +1, and Pfi = +1. Diamonds, 
representing 1/J3, indicate the lowest-energy state in the sector 
with K = 0, p,r = — 1, and pfi = — 1. Circles indicate other 
eigenstates in the equal-population case Nf = iVj. = N/2. 
Crosses show eigenstates for the minimally imbalanced case 
Nf = N/2+1 and N i = N/2— 1. Only the two lowest energies 
are displayed in each sector. 



factor, and then such factors are multiplied together. Fi- 
nally, we apply the symmetrization operation 6> groU p over 
all different ways of dividing the particles into k groups. 
Our motivation to consider the states (^) stems from 
the analogy with the scalar case; their single-component 
counterparts, Read-Reazayi [SU(2)k] states ||, give good 
approximations to the ground states of the scalar Bose 
gas H . The wave function (^) is a unique zero-energy 
eigenstate of a Hamiltonian consisting of a (k + l)-body 
interaction Eql: 



«i 2 )-"*(«<k-«ik+i) 5 ( 5 ) 



where the indices i±, . . . , ik+i run over all the particles 
of both components. On a torus, this Hamiltonian leads 
to a ground-state degeneracy of (fc + l)(fc + 2)/2 [ jl6| . 
The appearance of this non-trivial degeneracy is related 
to an underlying topological order, and can be used as 
an indicator of the corresponding topological phase. 

To discuss the possibility of SU(3)2 state at v = 
2/3 + 2/3, we present the energy spectra versus the pseu- 
domomentum K in Fig. ^j. At K = 0, we further decom- 
pose the Hilbcrt space using quantum numbers p w = ±1 
and pfi = ±1 for the 7r spatial rotation R n and the in- 
terchange of two components, Pfi, respectively. In ad- 
dition to the equal-population case Nf = JVj. = N/2, 
we also present data for the minimally imbalanced case 




0.5 



3 — r m _ _ ^ — a — e — A I 




— « — B^a — B — H — U — m — [ 

«/»]/,& SU(3) 2 □ ■ 




/ y 2 & SU(3) 2 o ~ 




(c) N v =6, N=8 ; 




0.5 0.6 0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 1 



L x /L„ 



U / U 



FIG. 3: (a,b) Energy spectra versus the aspect ratio L x /L y , 
at the filling factor v — 2/3 + 2/3 in the 5(7 (2)-symmetric 
case gfi — g. The same symbols as in Fig. are used. Solid 
curves are guides for the eyes for the states -01,2 • (c,d) Squared 
overlaps of -01,2 with the 517(3)2 states, plotted against the 
aspect ratio. 



find that for all the system sizes examined in Fig. the 
two lowest-energy states are spin singlets in the sector 
with K — 0, p n = +1 and pf± = +1 (labeled as ip\ and 
ip2)- By exact diagonalization of the 3-body Hamiltonian 
[Eq. (^|) with k = 2] on a torus, we find doubly-degenerate 
zero-energy eigenstates in the same sector. Therefore, 2- 
body and 3-body interactions seem to show a consistent 
topological degeneracy; due to a center-of-mass degen- 
eracy of 3, they both lead to the total degeneracy of 
6. Although the splitting between ipi and ip2 and the 
gap above them still show irregular dependences on N in 
Fig. ||, the consistency of the quantum numbers of these 
lowest-energy states with those of the 3-body Hamilto- 
nian provides an evidence for the formation of SU(3)2 
states. 

In Fig. ^ we examine the dependence on the aspect 
ratio L x /L y (only the case of L x /L y < 1 is examined 
without loss of generality). The ip\ and ip2 states con- 
tinue to be the two lowest-energy states in the ranges 
0.8 < L x /L y < 1 and 0.6 < L x /L y < 1 for N = 8 
and 12, respectively [Fig. ||(a,b)]. Such robustness of the 
spectral structures under a change in the ratio L x / L y is 
generally expected in vortex liquid states, which do not 
break the translational symmetry. We also calculated the 
overlaps of tpi and ip2 with the SU(3)2 states obtained 
as the ground states of Eq. (||) [Fig. ||(c,d)]. The average 
of the squared overlaps is around 0.8 and 0.6 for N = 8 
and 12, respectively in the range 0.7 < L x /L y < 1 (fo r 
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N/2 + 1 and N± = N/2 — 1 by cross symbols for L x /L y < 0.7, one of the overlaps decreases rapidly) 



N = 8 and 12; equal-population states which are not de- 
generate with crosses are identified as spin singlets. We 



Although the overlaps decrease as a function of TV, the 
increasing range of L x /L y with the robust spectral struc- 
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FIG. 4: (a) Energy spectrum versus coupling ratio g-fi/g 
for Nv — 6 and TV = 8 with L x /L y — 1. The states ^1,24 
and symbols are defined in a similar manner as in Fig. El 
The energy of ipi is subtracted from the whole spectrum, (b) 
Squared overlaps of the lowest-energy states Vi,2,3 with the 
states at g-f± = and the 5(7(3)2 states, plotted against g-fi/g- 



tures supports the appearance of the SU(3)2 state in the 
thermodynamic limit. 

Finally, we discuss a phase transition which occurs as 
the ratio g-fi/g is tuned at v = 2/3 + 2/3. For = 0, 
the system consists of decoupled scalar bosons; this sys- 
tem has been argued to show an incompressible com- 
posite fermion state at N a /Ny = 2/3 Q. On a torus, 
each component shows a triple center-of-mass degener- 
acy of ground states. Thus, the two-component system 
shows 9-fold degenerate ground states, each of which is 
given by the direct product of the above 3 states. In 
Fig. |](a), we present the energy spectrum as a function 
of g^i/g. For small g^i/g, we find three low-energy lev- 
els, which are separated from higher levels by a finite 
gap; each level is triply degenerate, implying the per- 
sistence of a quasi-degeneracy of 9 in this regime. For 



larger g^±/g (> 0.6), we find that the energy of ^3 goes 
up, while ipi and ip2 remain the lowest two states, lead- 
ing to a total quasi-degeneracy of 6 as expected for the 
5/7(3)2 state. These results point to a phase transition 
from the direct product of composite fermion states to 
the SU (3)2 state at g^/g ~ 0.6. Essentially the same 
behavior is found for a larger system size (Ny = 9 and 
N = 12; not shown), although the quasi-degeneracy of ipi 
and tjj2 around g\\,/g = 1 is less evident than in Fig. f|(a). 

The occurrence of the phase transition can be fur- 
ther supported by the results of the squared overlaps in 
Fig. ^(b). For small g^i/g, the states Vi,2,3 continue 
to have large overlaps with those in the decoupled case 
g^l = 0. For g^i/g(> 0.6), -01,2 have larger overlaps with 
the SU(3)2 states than with the states at = 0. It re- 
mains unclear whether a phase transition occurs directly 
between the two quantum Hall states or any intermediate 
phase exists. 

In summary, we have studied quantum Hall states in 
rapidly rotating two-component Bose gases. We have 
presented numerical evidences that a NASS state appears 
at v = 2/3 + 2/3. Non-Abelian anyonic excitations in the 
NASS state can carry spins (l5[ [l6| , and spin-selective 
operations would potentially offer better probe and con- 
trol of such excitations than in scalar Bose gases. Fur- 
thermore, we have demonstrated that changing the ratio 
g-fl/g drives a phase transition from the product of sepa- 
rate quantum Hall states to a spin-singlet quantum Hall 
state. Systematic appearance of spikes at v = fc/3 + fc/3 
with integer k in Fig. [j] suggests that SU(3)k states with 
k > 3 might also appear for the realistic contact interac- 
tions (|l|). More precise characterizations of the ground 
states at these filling factors deserve further studies. 
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Note added. During the preparation of this Rapid 
Communication, we became aware of an independent 
work by Grafi et al. p6| , where similar numerical evi- 
dences for the NASS states are presented. 
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